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Abstract 

We consider effects of a periodic modulation of the nonlinearity coef- 
ficient on fundamental and higher-order solitons in the one-dimensional 
NLS equation, which is an issue of direct interest to Bose-Einstein con- 
densates in the context of the Feshbach-resonance control, and fiber-optic 
telecommunications as concerns periodic compensation of the nonlinear- 
ity. We find from simulations, and explain by means of a straightforward 
analysis, that the response of a fundamental soliton to the weak pertur- 
bation is resonant, if the modulation frequency u) is close to the intrinsic 
frequency of the soliton. For higher-order n-solitons with n — 2 and 3, the 
response to an extremely weak perturbation is also resonant, if u is close 
to the corresponding intrinsic frequency. More importantly, a slightly 
stronger drive splits the 2- or 3-soliton, respectively, into a set of two 
or three moving fundamental solitons. The dependence of the threshold 
perturbation amplitude, necessary for the splitting, on ui has a resonant 
character too. Amplitudes and velocities of the emerging fundamental 
solitons are accurately predicted, using exact and approximate conserva- 
tion laws of the perturbed NLS equation. 

PACS numbers: 03.75.Lm, 05.45.Yv, 42.65.Tg 



1 Introduction 

The nonlinear Schrodinger (NLS) equation is a universal model of weakly non- 
linear dispersive media The existence and stability of solitons in the one- 
dimensional (ID) version of the NLS equation with constant coefficients is a well- 
established fact, which has important implications in various areas of physics. 
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In particular, solitons in fiber-optic telecommunications 3 and quasi- ID Bose- 
Einstein condensates (BECs) with attractive interactions between atoms 0], 
have drawn a great deal of interest. 

A new class of dynamical problems, which also have a vast potential for phys- 
ical applications, emerges in the investigation of soliton dynamics in extended 
versions of the NLS equation, in which coefficients are periodic functions of 
the evolutional variable. A well-known example is a nonlinear fiber-optic link 
subjected to dispersion management (DM), which implies that the dispersion 
coefficient periodically alternates between positive and negative values. The 
DM links support a family of stable temporal solitons (see, e.g., Refs. and 
also Ref. Somewhat similar is a waveguide- antiwaveguide system, which 

can be realized in the spatial domain (planar optical waveguides). In the latter 
case, a light beam is transmitted through a periodic concatenation of nonlinear 
waveguiding and antiwaveguiding segments 0. A common feature of the latter 
system with the DM is that the coefficient which periodically jumps between 
positive and negative values also belongs to the linear part of the equation. 

Another technique that may be useful for optical telecommunications is non- 
linearity management (NLM), which assumes that the coefficient in front of the 
nonlinear term periodically changes its sign. An advantage offered by the NLM 
is a possibility to compensate the nonlinear phase shift accumulating in pulses 
due to the Kerr nonlincarity of the optical fiber [Sj. In practical terms, the 
NLM can be implemented by dint of elements with a strong quadratic (x^) 
nonlincarity, which are periodically inserted into the fiber link. The \^ e l e ~ 
ments can emulate a negative Kerr effect through the cascading mechanism [Jj]. 
Various other schemes of NLM in fiber-optic links were considered, including 
its combination with the DM, amplifiers, etc. [U]|- A related scheme makes use 
of the NLM in soliton-generating lasers based on fiber rings [IT] . The NLM for 
spatial solitons, which assumes alternation of self-focusing and self-defocusing 
nonlinear layers in planar |12| or bulk |13j waveguides, was introduced too. 

All these systems may be regarded as examples of periodically inhomoge- 
neous optical waveguiding media. Other examples belonging to the same gen- 
eral class are tandem waveguides (see Ref. |14| and references therein) and the 
split-step model JS]. These are built as a juxtaposition of linear segments al- 
ternating with ones featuring, respectively, quadratic or cubic nonlinearity. A 
common feature of the models of all these types is that they support robust 
solitons, despite a "naive" expectation that solitons would quickly decay, pe- 
riodically hitting interfaces between strongly different elements of which the 
system is composed. 

The above-mentioned optical media are described by the NLS equation, in 
which the role of the evolutional variable belongs to the propagation distance, 
while the remaining free variable is either the local time (for temporal solitons) , 
or the transverse coordinate(s), in the spatial-domain models. Mathematically 
similar, but physically altogether different, models describe BECs in the ID 
geometry. In that case, the corresponding NLS equation is usually called the 
Gross-Pitaevskii (GP) equation. It governs the evolution of the mean-field wave 
function <f> in time (<), the other variable, x, being the coordinate along the 
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quasi-lD trap. In the normalized form, the GP equation is 

i<Pt = --0 KX + E/(a# + sl0| 2 </>, (1) 

where U (x) is the potential which confines the condensate, and the nonlinearity 
coefficient g is proportional to the scattering length of collisions between atoms. 
Two natural possibilities to introduce a time-periodic (ac) "management" in the 
BEC context are either through a periodic modulation of the confining potential, 
most typically in the form of U(x, t) = \ [kq + K\ cos(wi)] x 2 , or by means of 
time modulation of the scattering length, using the Feshbach resonance (FR) 
|16| . In the latter case, the nonlinearity coefficient in Eq. takes the form 
of g(t) = go + gisin(u>t). In either case, the modulation is generated by a 
combination of dc and ac magnetic fields applied to the BEC. 

The GP equation with the periodically modulated strength of the trapping 
potential was considered for both g > (when solitons do not exist, and the 
BEC as a whole is subjected to the "management", including the 2D and 3D 
cases) and g < 0, when the soliton is the basic dynamical object (TS|- In 
particular, a parametric resonance is possible in the former case, and creation 
of an effectively trapping potential, while the underlying one is anti-trapping, 
having n < 0, by the high-frequency ac part of the potential (with large ui) was 
predicted in the latter case. 

The periodic modulation of the nonlinearity coefficient through the "ac FR 
management" is an especially interesting possibility, as the FR is a highly effi- 
cient experimental tool, broadly used for the study of various dynamical prop- 
erties of the BECs ^Hj- In particular, it has been predicted that the modulation 
through the ac FR makes it possible to preclude collapse and generate stable 
soliton-like structures in 2D (but not 3D) condensates [T^l; m fact, this pre- 
diction is similar to the earlier considered possibility of the stabilization of 2D 
spatial optical solitons in a bulk waveguide subjected to the periodic NLM |13| . 
In the ID model of the GP type, subjected to the NLM, various stable dy- 
namical states, including Gaussian-shaped soliton-like objects, and ones of the 
Thomas-Fermi type, were studied in detail |2U) . In addition, analysis based on 
averaged equations was developed, for this case, in Ref. |2U (similarly to the 
analysis elaborated in Ref. ^H] for the case of the periodic modulation of the 
trapping potential). 

The objective of this work is to study resonance effects produced by the ac 
FR management, i.e., harmonic modulation of the nonlinearity coefficient, in 
the dynamics of fundamental and higher-order ID solitons in the NLS equation. 
We will focus on the case when the ac part of the nonlinear coefficient is small in 
comparison with its constant (dc) part go, which accounts for the self- attraction 
in the BEC, and is normalized to be g$ = —I. We also assume that the soliton's 
width is much smaller than the effective size of the trap, hence the external 
potential may be dropped. The consideration of the GP equation without the 
trapping potential makes it possible to identify fundamental dynamical effects 
for the solitons induced by the ac FR management. In this connection, it is 
necessary to mention that, in the ID case, the trapping potential is not a crucial 
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factor, on the contrary to the 2D and 3D cases, where the external potential 
plays a much more important role, in view of the intrinsic instability of the 
multi-dimensional NLS solitons. 

Thus, we will be dealing with the normalized NLS equation in the form [cf. 
Eq. Q] 

i*k + + [1 + fosin (cot)] |0| 2 = 0, (2) 

where the amplitude b of the ac drive is small. Note that Eq. (J2J) conserves 
exactly two dynamical invariants: the norm, which is proportional to the number 
of atoms in the BEC, 

/+oo 
\<f>(x)\ 2 dx, (3) 
-oo 

and the momentum, 

/+oo 
{Ml - 4>*<t> x )dx. (4) 
-oo 

We will demonstrate that the weak ac perturbation in Eq. (J2J can generate 
strong effects, if the driving frequency u) is close to specific resonant values. 
These effects include intrinsic vibrations of the fundamental soliton and splitting 
of the higher-order ones. We will also propose analytical explanations to these 
effects. To the best of our knowledge, these results have not been reported 
before for the present simple model. 

The rest of the paper is organized in the following way. The resonant effect 
of the periodic NLM on the fundamental soliton is reported in Section 2, and 
the resonant splitting of n-solitons with n = 2 and 3 is investigated in Section 
3. Section 4 concludes the paper. 



2 Resonant response of the fundamental soliton 

2.1 Numerical results 

First, we consider the action of the ac perturbation in Eq. (j2J on the funda- 
mental soliton, which, in the case b = 0, is 

<psoi(x, t) = A sech(A (x - x )) exp [iA 2 t/2) , (5) 

where A is an arbitrary amplitude. Numerical simulations of Eq. (0 ) were 
performed in a sufficiently large domain, < x < L, the initial condition cor- 
responding to the soliton (J5J placed at the center of the domain, xq = L/2. 
We have performed numerical simulations using the split-step Fourier method 
with 1024 Fourier modes. The system size is L = 50, and the time step for the 
numerical simulation is At = 0.001. 

Figure ^a) displays a typical example of the time evolution of the soliton's 
amplitude \4>(x = L/2)\, under the action of the ac perturbation with a very 
small amplitude, b = 0.0001. The frequency of the beatings observed in this 
figure can be clearly identified as u> — u> so \, where u> so \ = A 2 /2 is the intrinsic 
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frequency of the unperturbed soliton © . We have checked that the the beat- 
ing frequency is independent of the system's size and the other details of the 
numerical scheme. 
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Figure 1 : (a) A typical example of beatings in the evolution of the fundamental- 
soliton's amplitude under the action of a very weak perturbation in Eq. J3J, with 
b = 10~ 4 and uo — 2.2. The amplitude of the initial unperturbed soliton is A = 2 
(the corresponding soliton's frequency is lu so \ = A 2 /2 = 2). (b) The difference 
between the maximum and minimum values of the soliton's amplitude vs. the 
perturbation frequency. The dashed line is a fitting curve, 0.0003 • \u — 2|~ 1//2 . 




Figure 2: Direct numerical solution of the linearized equation © for the per- 
turbation around the fundamental soliton, in the near-resonance case, with 
oj = 1.98. 

The main resonant effect for the fundamental soliton is displayed in Fig. 
^b), in the form of the difference between the maximum and minimum (in 
time) values of its amplitude versus the driving frequency to. The resonance 
at u> = u! so \ = 2 (for A = 2) is obvious. The simulations do not reveal any 
noticeable subharmonic or higher-order resonance at frequencies u) = 1, 3 or 4. 
Due to the scaling invariance of Eq. (J2J , the plot shown in Fig. ^b) does not 
pertain solely to the particular value of the soliton's amplitude, A = 2, but is 
actually a universal one. It is easy to verify that the ranges of the variables t 
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and x, which are shown in this and other figures, correspond, in the normalized 
units, to experimentally realistic configurations of the BECs in the quasi-lD 
geometry. 



2.2 Perturbative analysis 

In order to explain the resonance shown above, we look for a perturbed funda- 
mental solution as <j>(x,t) = 4> ao \{x,t) + (j} pert (x,t), where the first term is the 
solution J5|- Thus, we arrive at the driven linearized equation for the pertur- 
bation, 

* (A»«t) t + \ W> P ert) ra + A 2 S ech 2 (Ax) (20 pcrt + e^'^crt) 

= l -bA 3 sech 3 (Ac)[e 4 (( Aa / a ) +w ) t -e < (( >lS / a )- w )*]. (6) 

The source of the resonant response is in the fact that the second term on the 
right-hand side of Eq. JSJ) becomes time-independent exactly at the resonance 
point, u> — uj so \ = A 2 /2. Figure|21displays the evolution of the perturbation close 
to the resonance at lu — 1.98, as found from direct numerical integration of the 
linearized equation with the initial condition 4>p e rt(x) = 0. As can be seen 
from the figure, the perturbation grows in time at the center, and simultaneously 
expands in space. Strictly speaking, the latter feature remains valid as long as 
the size of the region occupied by the expanding wave fields remains essentially 
smaller than the limit imposed by the confining field. 

The linearized equation 10 is too difficult for an exact analytical solution. 
However, the observation that the characteristic spatial scale of the solution 
observed in Fig. [21 becomes much larger than the internal scale of the function 
sech(^4a;) suggests that principal features of the solution can be understood 
from a simpler equation, in which the term ~ sech 2 (Ax) on the left-hand side 
of Eq. © is neglected, and the source corresponding to the second term on the 
right-hand side is approximated by a (^-function: 

i (</Wt) + 2 (^P Grt ) ra = COnst ' *( a! ) e " <Au ' t ) const = X 6 ^ 2 ' ^ 

where Aui = ui — ui so \, and const = (ib/2) A 3 sech 3 (Ax)dx. Equation 
can be solved by means of the Fourier transform. After straightforward manip- 
ulations, this yields 

~ . , const ,„ , A , , , f dt' ( ix 2 A ,\ 

<j> peI t(x,i) ^- — il + i^-^tj^ — expj^— +l A w -i'J . (8) 

Further consideration shows that, for Aw < 0, the asymptotic form of the solu- 
tion JSJ at t — > oo amounts to an exponentially localized stationary expression, 
which, by itself, is an exact solution to Eq. @: 

const 



f^pertOM) = ^==exp i|Aw|t- V2|A^N . (9) 
V 2 l Aw l v 
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In the case of Aw > 0, the asymptotic form of the general solution JSJ corre- 
sponds to a symmetric region occupied by plane waves emitted by the central 
source at wavenumbers k — ±V2AcJ. The region expands in time with the 
group velocities u gr = k = ±V2Aw, so that the asymptotic form of the solution 
is 



This asymptotic solution implies that the norm © of the expanding radiation 
field grows in time at the rate dN peTt /dt = y/2/Au> |const| 2 , which, in fact, is 
the rate at which the norm flows from the soliton to the radiation waves emitted 
under the action of the ac perturbation. 

Both analytical expressions 10 and (|10|l feature the lAwp 1 ^ 2 factor, that 
perfectly fits the numerical data summarized in Fig. ^b). Although these 
results, obtained for the weak time-periodic FR management, seem very simple, 
they have not been reported before, to the best of our knowledge. We also 
notice that the usual variational approximation (VA) for the NLS solitons, which 
is efficient in explaining a number of other perturbative effects [22], cannot 
account for the occurrence of the resonance at u> — uj so \, because the VA neglects 
radiation effects, while the above consideration showed that it is exactly the 
radiation field which is amenable for the manifestations of the resonance. 

The above results were obtained in the linear approximation, i.e., for a very 
small amplitude b of the ac drive in Eq. (J2J). At larger b, the perturbed soliton 
can either survive or decay into radiation. In fact, a stability region for the 
solitons in a similar model with a nonsmall perturbation, which differs from 
that in Eq. (J2J) by the form of the periodic modulation function, which is a 
piecewise-constant one, rather than harmonic, was drawn in Ref. |12| in the 
context of a model for spatial optical solitons in a layered waveguide. In the 
cases when a stable soliton established itself in the strongly perturbed ( "strongly 
nonlinearly-managed" ) system, its formation from the initial configuration JSJ 
went through emission of radiation and, sometimes, separation of a small sec- 
ondary pulse, while no pronounced resonance at u> = lu so \ was observed. As 
we do not expect that the replacement of the piecewise-constant modulation 
function by the harmonic one should dramatically alter the stability region, we 
do not consider this issue here in detail. 

3 Resonant splitting of higher-order solitons 

3.1 Response to a very weak ac drive 

As is well known, the unperturbed NLS equation gives rise to exact soliton 
solutions of order n, in the form of periodically oscillating breathers, which 
start from the initial conditions 




(10) 



4>o(x) 



nA sech(A (x - L/2)) 



(11) 
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with an integer n > 1 [the expression (|llf) assumes that the initial configu- 
ration is placed at the center of the integration domain] . The frequency of the 
shape oscillations (breathings) of the higher-order soliton is 

u; h , = AA 2 , (12) 

irrespective of the value of n. 

Generally speaking, the higher-order solitons are unstable bound complexes 
of fundamental solitons, as, in the absence of perturbations, their binding energy 
is exactly zero, which is a known consequence of the exact integrability of the 
unperturbed NLS equation. Nevertheless, not any perturbation readily splits 
the higher-order soliton into its fundamental constituents; usually, the splitting- 
is easily induced by specific nonconservative terms added to the NLS equation, 
such as the one accounting for the intra-pulse stimulated Raman scattering in 
optical fibers The consideration of dynamics of the higher-order solitons 
is also relevant, especially in the context of BECs, as the corresponding initial 
configurations can be created in the real experiment. 

We have studied in detail the n-solitons up to n = 5. First, we consider 
the case of a very small driving amplitude, b — 0.00005. Figure |31 (a) displays 
oscillations of the amplitude \(f>(x — L/2)\ of the 2 -soliton, which corresponds 
to the initial condition 1|11|) with n = 2. The frequency of the basic oscillations 
coincides with Wbr, as given by the expression (|12|l . while the frequency of the 
zoomed beatings in Fig. I^b) can be clearly identified with lo— uj^ x . The resonant 
character of the response of the 2-soliton to the weak NLM is obvious from Fig. 

13c). 



(a) (b) (0 




t t 



Figure 3: (a) Oscillations of the amplitude, \4>{x — L/2)\, of the 2-soliton created 
by the initial condition (flip with A = 1, in the case of b = 0.5 x 10~ 4 and 
to = 4.15. (b) Zoom of the previous panel around minimum values of the 
amplitude, which reveals beatings at the frequency lu — uj^. (c) The difference 
between the maximum and minimum values of the soliton's amplitude vs. the 
driving frequency uj. The dashed line is a fitting curve, 0.0006/|w — 4|; note the 
difference of the fitting power, —1, from that, —1/2, in Fig. ^b). 

Note that the exact solution for the n-soliton features not only the shape- 
oscillation frequency 1|12(1 . but also an overall frequency of the phase oscillations, 
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which coincides with the above-mentioned frequency o; so i = A 2 /2 for the fun- 
damental soliton, provided that the initial condition is taken as in Eq. 
In the simulations, we also observed a resonant response at lo = lo so \, but this 
resonance was essentially weaker than the one at to = Wbr- In particular, this 
is manifest in the fact that, as well as in the case of the fundamental soliton, 
the fit to the response around the former resonance is provided by the expres- 
sion \co — cj so i| , cf. Fig. H^b), while the fit to the resonance at lo = lo so \ 
demonstrates a more singular dependence, ~ \u> — u>t> r \ , as seen in Fig. Efc). 
Another qualitative difference between the two resonances is that the one at 
lo = Wbr, with a larger (but still small) forcing parameter 6, leads to splitting of 
the higher-order solitons into fundamental ones, as shown below, while, in the 
case of the resonance at lo = w so i, the increase of b does not lead to the splitting. 

3.2 Splitting of 2- and 3-solitons 

Unlike the case of the fundamental soliton, the reaction of the higher-order 
ones to larger values of the forcing parameter was not studied before, therefore 
we have done it here. First, we aim to demonstrate that the 2-soliton readily 
splits into two moving fundamental pulses, if the driving frequency is close to the 
resonant value (|12f) . The shape of each moving soliton is very close to that given 
by the commonly known exact solution, which can be obtained by application 
of the Galilean boost, with a velocity v, to the zero- velocity fundamental soliton 

cj) sol (x, t) = A sech(^(a; - vt)) exp [ivx + (i/2) (A 2 - v 2 ) t] . (13) 

Figure 01 displays the evolution of the wave function for the initial condition 
with n = 2 in the resonant case (^4 = 1 and lo = 4), with the driving 
amplitude b = 0.0005. The latter value is still very small, but larger by a factor 
of 10 than in the case shown in Fig. The amplitudes of the two fundamental 
solitons, observed as a result of the splitting, are close to A\ = 3 and A2 = 1 
[note that they exactly corresponds to the fundamental-soliton constituents of 
the original 2-soliton with A = 1, in terms of the inverse scattering transform 
(1ST) 221]. Velocities of the splinters were measured to be v\ = 0.00197 and 
i>2 = 0.0066, respectively (with the ratio V\ : V3 w 1 : 3). At the end of the 
simulation run (t = 1000 ), the secondary solitons are found at the distance, 
respectively, 4.5 and 13.2 from the central point, x = L/2. 

Similar near-resonant splittings were observed for n-solitons with n > 2. 
In particular, Fig. shows this outcome for n — 3, which corresponds to the 
initial configuration (|11|) with n = 3, A = 0.5, lo = 1 and b = 0.0005. This time, 
the splitting gives rise to three moving fundamental solitons, whose amplitudes 
are close to Ai = 2.5, A2 — 1.5, and ^3 = 0.5. As well as in the case of 
n = 2, these values correspond to the constituents of the original 3-soliton (with 
A = 0.5), in terms of the 1ST [23] - The velocities of the three splinters are 
vi = —0.00146, v 2 = 0.0732, and ^3 = —0.0148, so that the ratios between them 
are V1/V2 ~ —1/5 and V3/V2 ~ —2. 
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Figure 4: A typical example of the splitting of a 2-soliton [generated by the 
initial condition with n — 2 and A = 1] into an asymmetric pair of moving 
fundamental solitons, under the action of the weak resonant drive, with uj = 4 
and b = 0.0005. (a) The evolution of |w(a;,t)|. (b) The final configuration at 
t = 1000. 




Figure 5: The same as in Fig. 01 for the 3-soliton, generated by the initial 
configuration (|llfl with n = 3 and A — 0.5. In this case, the forcing frequency 
and amplitude are ui = 1 and b = 0.0005. 
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These results can be summarized in the form of diagrams which show the 
minimum (threshold) value of the forcing amplitude b, necessary for the split- 
ting, versus the driving frequency u>. The splitting of the 2- and 3-solitons was 
registered if it took place in the simulations of Eq. J2J that were run up to the 
time, respectively, t = 600 or t — 2000 (still longer simulations did not give 
rise to any essential difference in the results). As is seen from Fig. for both 
2- and 3-solitons these dependences clearly have a resonant shape, with sharp 
minima at the frequency given by Eq. (|12|l . It is not quite clear why the forcing 
amplitude required for the splitting is very small but finite even exactly at the 
resonance point. This may be related to the accuracy of the numerical scheme 
and/or the finite size of the integration domain. Similar observations were also 
made in simulations of the n-solitons with n — 4 and 5. 

(a) (b) 

0.003 | 1 1 1 1 1 — 

0.0025 k 
0.002 - \ 

0.0015 - \ y- 

0.001 - \ y 
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I ' ' ' ' ' — 

3.8 3.9 4 4.1 4.2 0.9 0.95 1 1.05 1.1 




Figure 6: The minimum values of the amplitude of the ac perturbation, neces- 
sary for the splitting of the 2-soliton (a) and 3-soliton (b), as functions of the 
driving frequency. The initial condition is taken in the form of Eq. (|1 1|) with, 
respectively, n — 2 and A = 1, or n = 3 and A — 0.5. In both cases, the sharp 
minimum exactly corresponds to the resonant frequency, as predicted by Eq. 



3.3 Analytical results 

The amplitudes and velocities of the fundamental solitons, into which the higher- 
order ones split, can be predicted in an analytical from. As it was already 
mentioned above, the amplitudes of the secondary solitons coincide with those 
which correspond to the constituents (eigenvalues) of the corresponding original 
n-soliton in terms of the 1ST. However, the velocities of the emerging funda- 
mental solitons cannot be forecast this way, as, in terms of the 1ST, they are 
zero when the fundamental solitons are bound into a higher-order one. 

Nevertheless, both the amplitudes and velocities of the final set of the solitons 
can be predicted in a different way, using the exact and nearly exact conservation 
laws of Eq. @. Indeed, there are two exact dynamical invariants, J3J) and (QJ, 
and, in addition to that, the unperturbed NLS equation has an infinite series of 
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higher-order dynamical invariants, starting from the Hamiltonian, 




(14) 



Two next invariants, which do not have a straightforward physical interpreta- 
tion, are |2] 



In the case of the splitting of the 2-soliton Hll|) with the amplitude A, the 
exact conservation of the norm and approximate conservation of the Hamilto- 
nian l|14l) yield the following relations between A and the amplitudes Ai$ of the 
emerging fundamental solitons (splinters): 4A = A\ + A2, and 28 A 3 w A\ + A\ 
(the latter relation neglects small kinetic energy of the emerging solitons) . These 
two relations immediately yield A\ = 3A and A2 — A, which coincides with the 
the above-mentioned numerical results, as well as with the predictions based 
on the set of the 2-soliton's 1ST eigenvalues. Furthermore, the exact momen- 
tum conservation yields a relation involving the velocities of the secondary 
solitons, A\V\ + A2V2 = 0. With regard to the ratio Ai/A 2 = 3, this implies 
v\/v2 — —A2/A1 = —1/3. This relation is indeed consistent with the aforemen- 
tioned numerical results, although the absolute values of the velocities cannot 
be predicted this way. 

Similarly, in the case of the splitting of the 3-soliton, the exact conservation 
of N and approximate conservation of H and I5 [see Eq. ()15|l ] yield the relations 
(which again neglect small kinetic terms, in view of the smallness of the observed 
velocities) 9A = Ai+A 2 +A 3 , 153A 3 w A\+Al+A% and3369A 5 w Af+A|+A|. 
A solution to this system of algebraic equations is A\ = 5A, A2 — 3A, A3 = A, 
which are the same values that were found from the direct simulations, and can 
be predicted as the 1ST eigenvalues. The conservation of P and I4 gives rise to 
further relations, Ai^i + ^2^2 + ^3^3 = and {A\v\ - A\ v\) + (A 2 w| - A^) + 
(Aawl — A3W3) = 0. If the velocities are small, it follows from here that 



Vl /v 2 = -{Al-A 2 Al)/{A\-A l AD = -1/5, and v 3 /v 2 = -{Al-A 2 A\)/{Al- 



A3AI) = —2. These results for the velocities are consistent with the numerical 
situation observed in Fig. 5. 

4 Conclusion 

In this work, we have addressed a simple model, based on the NLS equation, 
which describes an attractive Bose-Einstein condensate (BEC) in a quasi-lD 
trap, with the nonlinearity strength subjected to a weak time-periodic (ac) 
modulation (that can be imposed by means of the Feshbach-resonance tech- 
nique). The same model describes the nonlinearity management in periodically 
inhomogeneous optical waveguides. 





(15) 
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It was found from direct simulations, and explained by means of a straight- 
forward perturbative expansion, that the response of a fundamental soliton, 
in the form of temporal beatings of its amplitude, to the weak ac perturba- 
tion is resonant when the driving frequency u is close to the soliton's intrinsic 
frequency. For n-solitons (breathers), with n — 2 and 3, the response to an 
extremely weak drive is also resonant, if w is close to the breathing frequency. 
More interestingly, a slightly stronger drive gives rise to splitting of the 2- and 
3-solitons into sets of two or three moving fundamental solitons, respectively. 
The dependence of the minimum perturbation amplitude, which is necessary 
for the splitting, on u) has a clearly resonant character too. The amplitudes of 
the splinter solitons, and the ratio of their velocities, can be easily predicted on 
the basis of the exact and approximate conservation laws of the perturbed NLS 
equation. 
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